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EQUATIONS WITH (1 + 1) POTENTIALS 
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Abstract. We study the Darboux transformation (DT) for Dirac equations with (1 + 1) poten- 
tials. Exact solutions for the adiabatic external field are constructed. The connection between 
l/^ , the exactly soluble Dirac (1+1) potentials and the soliton solutions of the Davey-Stewartson 

' equations is discussed. 

o ■ 

(N 



O 
(N 



> 



1. INTRODUCTION. 

The Darboux transformation (DT) is a convenien way to construct a rich set of integrable 
potentials of the steady-state Shrodinger equation in the single dimensional case [1]. As 
shown in [2], DT can be used to study the Dirac equation for a two-dimensional fermion 

■ in an external scalar field w{x). The aim of this work is the generalization of DT for the 
0> , Dirac equations with one-space-dimensional and non-stationary potentials u{t,x). In Sec. 
^ ! 2 we show that for a fermion of transverse momenta p = py, q = the four-dimensional 
O I equation reduces to the Zakharov-Shabat equation. We demonstrate DT for this equation 
^ ■ and the connection with intertwining and supersymmetry algebra. The result of the multiple 

■ DT (extended Crum law [3]) is set forth in Sec. 3. In Sec. 4 we show the connection between 
the exactly solvable Dirac (1+1) potentials and the soliton solutions of the Davey-Stewartson 

Qh! equations [4]. We also discuss the reduction restriction problem. 
^ ; 2. DARBOUX TRANSFORMATION 

Let us consider the four-dimensional Dirac equation 

ih^d^ - rMt^ x) + m)* = 0. (1) 
We use the 7-matrix representation [2] , [5] 
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and for A'^{t, x) we have either: 

A^it, x) = (0, 0, A^{t, x),A^{t, x)f, (3) 

or 

A''(t,x) = (0,0,(5(t,x),cQ(t,x))^, c = const. (4) 
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One can easily see that (3) and (4) are solutions of Maxwell equations. Let ^ have the form: 

^ = ^{t,x) exp{i{py + qz)), Imp = Imq = 0. (5) 



Then $ satisfies the Zakharov-Shabat equation: 
For the case (3): 



U = 



where A = i{A^ - q), B ^ i{p - A'^). 
For the case (4) we get: 



/ 





A 





m — 


B\ 




A 


^ 


B-m 












B + m 





A 




V 


-B-m 





A 





J 



/I \ 
0-100 

10 

\o -1/ 



where 
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iq{n^ + if) + 2mn3 
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The condition c* = c gives us for a = + iaj , /? = + iP^ 

+2maj{(3^'^ + a^'^){pa^ + q(5j^) = 0. 
One of the nontrivial solutions of (10) is 



i{ap + /%) + am 



(3, =0, (3,= 



a 
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(6) 
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Substituting (8), (9) into (1) and taking (10) into account, we get a 2 x 2 equation (6) for 
F where J = (73, 

'0 u 
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u{t, x) = \i + MQiti x), v{t, x) = Vi + V2Q{t, x), 



Ai = 



m/3 + i{qa — p(3) 
1^ 



1 (a^ + P'^){qa - p(3 - im(3) 



_ ^{rrf + q"^ + p"^) _ '^'"n ^ 

i{p/3 — qa) — m(3 ' (3 



(12) 



Let Xi and X2 be 4 X 4 (for the case (3)) or 2 x 2 (for the case (4)) matrix solutfons of 
equation (6). We define a matrix function = Xi,xX^^- It ^asy to see that satisfies the 
foUowing nonhnear equations: 

^i,t = ^sr,,, + [U, rj + [as, + U^. (13) 

Equation (6) is covariant with respect to DT: 

X.[l] = X2,x-r,x„ U[l] = U+[a3,T,]. (14) 

It is necessary to choose the function Xi in such a way that the structure of the matrix 
U[l] be the same as the structure of the matrix U. This is the condition that we call the 
reduction restriction (see Sec. 4). 

The transformation (14) allows us to constract a superalgebra, in just the same way as the 
DT for the steady-state Shrodinger equation in a one-dimensional case [6]. In order to do 
this we introduce the following operators: 

Let us define new several operators as follow: 

h = G^-^G^+\ h[l]=G^+^G^-\ (16) 

It easy to see that 

G(+^T = T[l]G^+\ TG^-) = G^-)r[l], [h,T] = [h[l],T[l]]=0. (18) 
The operators h and h[l] are the tipical one-dimensional matrix Hamiltonians: 

P^=(rt-r,)u, V = -(t,^,+t+t,), ^[1] = r+ - r,r+, (20) 
where (rj^ - diagonal part of r^. It easy to see that the operators q^^\ H 

«=(j.> °), (21) 



« -\G(+i Oj' « -1,0 J' "-{O h[l] 

generate a supersymmetry algebra: 

= [q^^\H] = 0, {q^+\q^-^} = H. (22) 

Furthermore we restrict ourselves to studing the case (4). The case (3) will be considered 
in a separate work. 

3. EXTENDED CRUM LAW 
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Let us consider 2N+1 particular solutions of (6) with = (ipk, 0^) , k < 2N, $ = (ip, 4>): 

i'k,t + iJk,x+u{t,x)4>k = 0, (l)k,t- (f>k,x + v{t,x)ilJk = 0. (23) 

The following theorem is established: 
THEOREM 

Functions ip[N], (f)[N] satisfy (23) with potentials u[N] and v[N] such that: 
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u[N] = w + 2 
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(24) 



where Ai^2, -Di,2, D are the following determinats 
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AA: 



02JV'' V'S ^) - V'2 
To prove this theorem we construct N 2x2 matrix functitons Xk — (^2fc-i, ^2fc), 1 < ^ < 



^•2fe-l ^32fe 
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These functions satisfy the matrix equation (6) with J = as- 
After N-time DT (14) we get x[N] and U[N] which satisfy (6). Let us write xW] as a 
series: 



X[N] = x(^) -J2Mt,x)x^^-% A, = h M 

i—l ^ ' 

Plugging (25) into (23) we get (here = kiiN-ky. )'- 



(25) 



AT N 

k=0 fc=l 

AT N-k N 

+ E E Cir^feC/^X^^-'^-^^ - U[N] (x(^) + E ^feX^^-'^) = 0, (26) 

fe=l i=0 fe=l 

where (^fe)^- is the off-diagonal part of A^. Therefore 

u[N]^u + 2bi, ^;[A^] =i;-2ci. (27) 

To compute bi and ci we take into account that XA;[-^] = Q if k < N, therefore we get a 
system of 2N equations as follows: 

N N 
n=l n=l 

i — 1, N. Using Kramer's formulae, we get (24). 

If (see(12)) |U = — ca*)~^ then Imv2 = lTn\2 = 0. The constants Ai and vi may 
be annihilated by the standard gauge U{1) transformations. Using freedom in our choice of 
parameters, let us assume that v(t, x) = Ku{t, x), k = ±1. Now we can supplement (24) with 
the transformations law for Xi {N = 1): 

I rn ^2 I rn rn ^^^^ 



^1 ^2 

01 02 



(29) 



It is necessary to require that after N-time DT the following reduction restriction will be 
true: 

v[N] = Ku[N]. (30) 

In the general case we do not have an algorithm allowing us to keep (30) in all the steps of 
DT. However, it is possible by the introduction of the so called binary DT which allows one 
to preserve the reduction restriction (30) [7] . 
Let us consider a closed 1-form 

*. = <ixc, + *ca3x, n.Jin (31) 

where a 2 x 2 matrix function ( solves the equation: 

Ct = - CU. (32) 
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We shall apply the DT for (6). One can verify by immediate substitution that (32) is covariant 
with respect to the transform if 

C[+1] = ^C,X)X-'- (33) 
Now we can alternatively affect (7, by the following transformaton: 

U[+l,-l]^U + [as,xn-\]. (34) 

It may be shown that 

TV N 

X[+N, -N]=x-Yl ^fe^(Cfe, X), C[+N, -N]=C-Y1 ^(C> Xk)sk, (35) 

fc=l k=l 

where $k and Sk may be found from the following equations: 

N N 

^Ok^{Ck,Xi) ^ Xi, ^^{Ci,Xk)sk ^ Ci- (36) 

/c=l k=l 

The transformation: 

N 

U[+N, -N] = U+Y, [^3, OMCk, Xi)sk] (37) 



i,fe=l 

is the forementioned binary DT. 

Let V — Ku* (see Sec. 4; in this section u and v are real so the condition is equivalent to 
(30)), then U[+N, —N] will satisfy the reduction restriction if we choose (k and Xk such that: 

Cfe = XkR, R = diag{l, -k). (38) 
The solution that follows from (34) has the form 

u[+l,-l]=u+ — , (39) 

Oik ^ J dx {ipiipk - ii4>*i4>k) + dt {ip*ipk + n(l)i4>k)- 
Note that the square of the absolute value u[+l, —1] is expressed by the compact formula: 

I u[+l, -1] |2=| u f -k{^ - ^) ln(^n^22- I 6,2 \'). (40) 

DT allow us to constract a rich set of the exact solutions of the Dirac equations with (1+1) 
potentials. In conclusion of this section we consider decreasing in t, x — > ±oo fields (adiabatic 
engaging and turning-off). Let u = v = 0. The particular solutions of (23) are: 

ijk = Ak e'^^*--) + Bk e'^(--*), 0^ = Ck e^^*+") + Dk e-^^'+''\ (41) 

where k = 1,2; A, B, C, D, to and A are real constants. After DT (14) under condition (12) 
we get: 

Q(t^x) = - , (42) 

' Hi cosh(ai + bx + Si) + fj,2 cosh(6t + ax + 62) ' 
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1 + 

u{t,x) ^ ^Q{t,x), v{t,x) = ■^u{t,x), (43) 

^ =1 ac - /3 p== 4u;{AiB2 - A2B1), a = u + \, b = X-u, )Lt(/3* - ca*) = i, 
1 + c2 = 4Ae(C2L'i - CiDa), uj{AiB2 - A2S1) = X{DiC2 - D2C1), 



f^i = V{AiC2 - A2Ci){BiD2 - B2D1), 1x2 = V{BiC2 - B2Ci){AiD2 - A2Di), 

oi — - In , Oo = - 111— : , 

2 B1D2-B2D1' 2 AiD2-A2Di' 

where a, P satisfy (11). The potential (42) is a locahzed impulse, which decreases in t — > 
±00. A particular solution of the Dirac equation may be calculated using (5), (8) from 
('^i,2[l],0i.2[l])"^, where V'i,2[l] and (/>i,2[l] are defined by (29). This bispinor describes the 
quazi-stationary state of a fermion. The fermion is free along the y and z and restrained 
along the x-axis. 

4. DAVEY-STEWARTSON EQUATIONS. 

It is shown in [2] that a class of the exactly soluble Dirac one-dimensional potentials 
corresponds to a soliton solutions of the MKdV equation, just as certain Schrodinger potentials 
are solutions of the KdV equation. In this section we show that there exist the analogous 
connection between exactly soluble Dirac (1+1) potentials and soliton solutions of the Davey- 
Stewartson equations (DS). 

The DS equations appear as the commutation condition of the two operators [7] , [8] : 

[Ti,T2]=0, (44) 
where Ti = T (see (17)) and T2 is given by the formula: 

T2 = i^+ 2as^ + 2U— + U, + a^Ut + A, (45) 
ay ox^ ox 

A = diag{Ai, A2), v = ku* , 

1 1 
Ai = -K\uf +Yik+)^^ A2 = k\u\^ +2i^^")^' 

Here /(j-) — and -F is a pure imaginary function. The operator T2 is also covariant with 

respect to DT, what allows one to get infinite sets of exact solutions of DS, for example soliton 
solutions, exponentially localized on the plane - the dromions (42) [7], [8]. The DS equations 
contain two fields: u{y, t, x) - the amplitude wavetrain function and S{y, t, x) = —iFx - the 
amplitude function of slowly changing in t (space variable!) and x mean field: 

iUy + Uxx + Utt — 2K\u\'^U+Su=0, = -4k(| It |^)a;a:. (46) 

The dromions iS-component is constant along the two ortogonal directions and moves with 
the constant velocity [8]. 

In Sec. 3 we have shown that the one-dromion solution corresponds to adiabatic engag- 
ing and turning-off electro-magnetic field. Boity, Leon, Martina and Pempinelli [8] showed 
that two dromions scatter against each other with the soliton phase shift, therefore the N- 
dromions solution corresponds to a adiabatic external field, too. A many-dromion solution 
may be obtained by the N-time DT. For this purpose we transform the LA-pair (Ti and T2). 
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Introducing new variables x = p + q and t = p — qwe exclude a field F from (45) . The matrix 
A takes the form 

where gi{p,y) and g2{q-iy) are arbitrary functions. We suppose that u{p,q,y) belongs to a 
Schwarz space L. Then a nonlocal flow is given as follows: 

5 = 2k I « |2 +( dq^+ r dp-^)\uf+gi- 92. (48) 
J-00 op J_oo aq 

In the new variables ^ 

therefore for u & L, the equation Tix = gives at infinity for components of x the following 
condition: 

fpk = ipk{p,y), (t>k = (t>k{q,y), k = i,2. (50) 

The second equation T2X — gives four nonstationary one-dimensional Schrodinger equations 

#fc,y + \'^k,pp + Viljjk = (51) 



#fe,2/ - ^^k,qq + V24>k = (52) 



where 



MP^y) ^ 91 + J dQ : V2{q,y) ^ 92 + 1^ J dp . (53) 

For the case of a dromion (42) we have 

vi = Auj^sec^fi, V2 = -4A^sec^/2, 

/i = 2u{p - 2ay) + ^ In h= 2A(g + 2hy) + ^ In ^, (54) 

2 bo 2 Do 

where a, 6, {A, B, C, D)o are real constants. So the dromion sits on two plane solitons that 
correspond to the one-level refiectionless potentials. They may be easily obtained by the 
standard DT for the Schrodinger equation [1] on zero background: 



, rn dp '^0 gp iPk Qq (PO Qq (Pk 



^i[l] = -2^1nV^o, ^2[l] = 2^1n0o. (55) 



The background eigenfunctions are choosen such that: 



V'o (P, 2/) = 2 ^/AoBo cosh /i e'^ , (56) 
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Mq,y) = 2v/Co^cosh/2e^^ (57) 

where 

d = 2[{lu^ - a'^)y + ap], ^ = 2[{b^ - X^)y + bq], (58) 

and the new coefficients (A, B, C, D)o satisfy the conditions: 

AoiB^C^ - B^C^) = Bo{A^C2 - A^C^), CoiA^D^ - A^D^) = DoiB.C^ - B^C^). 

For the second DT it is necessary to switch to the two-level potentials in equations (51), 
(52). Let w > uj, g > X then two "linearly independent" solutions that correspond to the 
eigenvalues 2zu and 2q have the form: 

for (51) and 

= C_i e2^'('?+2/5y)+i^-i ^ 0(_-) ^ e-2^('^+2^^)+^''-i (60) 

for (52), where 9-i and -d-i differ from 6 and by substitution (a, fe. A, uj) {a, (3, q, uj). We 
transform and ^^'^ by the formula (55) with ipo and 0o accordingly. Then we will build 
the new support function: 

v^_i[i] = - [1], <^_i[i] = 0L+i)[i] - 0y [1]. (61) 

For simplicity let us choose A = B = C = D = 1. Asa result we have: 

^_,[1] = ( - ^) cosh(^o + Ci) + (^ + ^) cosh(eo - Ci) ^ .^^ _ ^) ^.^1^ ^g2) 

cosh 4o 

</>_,[l] = ( (g-A)cosh(r;o + ryO + (, + A)cosh(r,o-^0 ^ .^^ _ ^.^^^ ^^^^ 

cosh ryo 

where = 2u!{p — 2ay), = 2tu(p — 2ay),r]o — 2X{q + 2by), rji — 2g{q + 2by). If one puts 
a = a, P = b then (62), (63) exactly coincide with the supporting functions that generate the 
two-level potential with respect to the Darboux transformations [8] . 

Now it is necessary to define four functions V'fcfl], 0fc[l] {k — 1,2) which are the solutions 
of (51), (52) with the potentials (54). The basic problem at this step is the independence of 
the wronskians 
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with respect to the space coordinates, what is necessary for the solvability of the reduction 
restriction equation v[2] = ku*[2]. For this purpose it is enough to choose: 

tanh ^o)ip^_X^ + Bk-iiw + u tanh ^o)V'ii\ (64) 

<Pk[l] = Ck,-i{g - X tanh Vq)(P^-i + Dk,-i{g + X tanh f]o)(p^~i\ (65) 

with the constants that satisfy the relation: 

w{cv^ - w^){Bi,.iA2,-i - 52,-iAi,_i) = KgiX^ - g^){C^,-^D2,-i - C2,-iZ?i,-i). (66) 

Therefore we choose ?/;_2[2], 0-2 [2] as the wave functions that generate a three-level potential 
(via DT) and determine the corresponding functions '0fc[2] and 0fe[2]- One may repeat this 
procedure and realize the third DT. 
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